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In the present paper we have analysed a fermionic infinite-ranged quantum Heisenberg spin glass 
(s = 1 /2) with a BCS coupling in real space in the presence of an applied magnetic field. This model 
has been obtained by tracing out the conducting fermions in a superconducting alloy. The magnetic 
field is applied in the resulting effective model. The problem is formulated in the path integral 
formalism where the spin variables are defined as bilinear combinations of the Grassmann fields. 
The static approximation is used to treat both the pairing and the spin terms together with the 
replica symmetry ansatz. Henceforth, the problem can be reduced to a one site problem. The field 
in the z direction, H z , separates the order parameters in two groups: parallel and transverse to it. 
We have obtained a phase diagram in T—g space with zero transverse spin-glass ordering, g being the 
strength of the pairing interaction. It has been possible to locate the transition temperature between 
the normal paramagnetic phase (NP) and the phase where there is a long range order corresponding 
to formation of pairs (PAIR). The transition ends at the temperature Tf, the transition temperature 
between the NP phase and the spin glass (SG) phase. Tf decreases for stronger fields allowing us to 
calculate the NP-PAIR line transition even at low temperatures. The NP-PAIR transition line has 
a complex dependence with g and H z , having a tricritical point depending on H z from where second 
order transitions occur for higher values of g and first order transitions occur for lower values of g. 
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I. INTRODUCTION 

The interplay between a magnetic order and superconductivity is a current issue. Particularly, spixglass ordering has 
been reported in many physical systems that includes high-T c superconductorscm, heavy fermionsfju and conventional 
superconductors doped with magnetic impuritiesO like Gd x Th\_ x RUi. These conventional superconductors can be 
modeled by a s-d exchange interaction betweeruthe localized magnetic impurities and the conducting fermions which 
have conventional BCS interaction (Nass et alfl). The main interest of Nass et al. was to find the density of states 
in the presence of magnetic impurities. However, relativeLyJittle consideration has been given to the phase transition 
problem between spin glass and a superconductive orderuEl. 

Nevertheless, from a wider perspective the role of quantum fluctuations in disordered strongly interacting systems 
still needs to be fully understood. Theoretical studies in recent years have shown plenty of interesting results. For 
instance, the presence of a non-Fermi liquid behavior near T = transition between a metallic paramagnetic and-a 
metallic spin glass in a model introduced to study the competition betweea-Kondo effect and RKKY interactiontll. 
In addition, the random quantum Heisenberg with generalized SU(M) spinal^l produces a spin liquid ground state in 
the large M limit with a local dynamic susceptibility that is identical to a marginal Fermi Liquid. 

Recently, a model has been introducecLtp study the phase transition between the spin glass ordering and the BCS 
pairing among fermions of opposite spinsE3. This model has been obtained from the model of Nass et al.El by tracing 
out the conducting fermionic degrees of freedom by perturbation expansion. The remaining effective problem could 
be solved in its infinite ranged Ising version. A rich phase diagram has been found in the T—g plane (g is the strength 
of the pairing interaction) for fixed J (variance of the spin random coupling «7y) with a normal paramagnetic phase 
(NP), a spin glass phase (SG) and a pairing phase (PAIR) where there is formation of pairs in the sites. The first two 
phases (NP and SG) are separated by a second order line at the spin glass transition temperature Tf. However, the 
line g = g c (T) between the NP and PAIR phases is more complex. It is a second order line which ends at a tricritical 
point (T tc ,g tc ), where a first order line starts. Below the spin glass transition temperature Tf, the last two phases 
(SG and PAIR) are separated by the same first order line transition. 
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The purpose of the present study is to investigate the infinite ranged Heisenberg version of the model given by 
Ref. [12| in presence of a magnetic field for the half filling case. This many fermion problem is formulated by defining 
the SU(2) spin variables as a bilinear combination of Grassmann fields and computing the partition function through 
the path integral formalism. The disorder in the coupling is treated by standard replica formalism. For this 
representation of the spin variable, the Fock space (s — 1/2) has four quantum states per site but only two of them 
are magnetic. As a consequence, the sensitivity of the site to magnetic coupling is dependent on the occupation 
number which is controlled by quantum statistics. When the strength of the BCS coupling is increased, forcing the 
double occupation, the number of sites non-sensitive to the magnetic coupling increases. Therefore, this mechanism 
could suppress magnetic ordering such as the spin glass phase. Actually, that is what happens in the Ising modelij 
even at the half filling situation. 

In the replica formalism for the quantum spin glass, the replica diagonal order parameter has to be solved coupled 
with the non-diagonal one (the spin glass order parameter) and the remaining order parameters of the problem. In 
consequence, the diagonal component is no longer constrained to unity and has an important role in the problem. In 
fact, the location of Tf can be obtained through this component. This is the essential point: the location of other 
transitions of the problem are also dependent of this diagonal component even when the non-diagonal component is 
zero. 

Another important aspect of these problems is that the order parameter (the spin-spin correlation function) are 
time dependent. Bray and More£3, in their pioneering work, investigated the role of the quantum fluctuations in 
the Heisenberg spin glass model by using the static approximation. They formulated the problem in the Feynman's 
path integral formalism where the spins variables have a dependence on fictitious time r. The static approximation 
neglects the time dependence giving an upper bound for the thermodynamic potential and the location of the transition 
points. For temperatures above the Tf , this approximation is expected to be exact together with the replica symmetric 
solution. This same problem has been studiedt-3 with a formalism where the spin variables were represented by bilinear 
combinations of Grassmann fields. The static approximation was used for the spin glass order parameter, while for 
the susceptibility (associated to the replica diagonal order parameter) was applied an instantaneous approximation. 
This approach has found a non-zero spin glass transition temperature for s = 1/2. Goldschmidt and Pik-Yin Lali3 
used the same treatment as Bray and Moret3 to obtain a phase diagram for the infinite ranged quantum Heinsenberg 
and XY spin glass in an external magnetic field. Due to the magnetic field, the order parameters were separated 
in two groups: transversal to the field and parallel to it. This approach allowed to locate the spin glass transition 
transversal to the field and investigate the role of the tunneling in the transition. For instance, the results showed 
that the spin glass (s = 1/2) transition temperature Tf is depressed for the Heisenberg model as the field strength is 
increased. Nevertheless, the tunneling was not strong enough to produce a transition at T = 0. 

In the present fermionic formulation, we have combined the methods of Refs. [l4| and [l5| by applying a field over the 
fermionic spin variables, taking the transversal spin glass order parameter as null and using the static approximation 
for the remaining order parameters. In these circumstances, the magnetic field also has depressed the transversal spin 
glass transition temperature. This approach has allowed us to study the transition line between NP phase and the 
PAIR phase (we use here the same terminology of Ref. ^) even in lower temperature. This transition line ends at 
the transversal spin glass transition temperature which moves down as the field strength enhances. 

This paper is organized as follows. In section [n[ we present our model and develop it to obtain the Gran Canonical 



Potential and the saddle point equations for the order parameters in the half filling case. In section IB, we analyse 



the transition line between the NP and PAIR phases for several values of the magnetic field. The results show that 
as the magnetic field increases, larger values of the pairing interaction strength are necessary to produce a transition. 
We also locate the transition line between the normal paramagnetic phase and the spin glass phase. Discussions and 
concluding remarks are presented in the last section. 



II. GENERAL FORMULATION 

We consider the following HamiltonianEl in the presence of a magnetic field H z 



i s=T,i i,3 i 



where af = af ± of, the chemical potential is given by \i and cj s (cj s ) are fermions creation (destruction) operators 
(s =| or J, means the spin projection). The spin operators are represented as bilinear combination of fermionic 
operators: 
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u t = \\ c \\ c a - c !i c a] (2) 

(jf = 4 ]Cii . (3) 

The random coupling Jij in Eq. (Q) are infinite ranged with a Gaussian distribution with mean and variance given 
by: 

(4> ca = AJ 2 /N (4) 

where N is the number of sites. 

The third term in Eq. (Q) represents a BCS pairing interaction in. r e al space. This pairing mechanism in fermionic 
problems already has been purposed, although, in different context! 1 3il3. The remaining terms describe a fermionic 
Heisenberg spin glass in the presence of a field. Our interest in this fermioniC|-Broblem is to analyse the phase 
transitions at mean field level as firstly introduced by Sherrington and Kirkpatricku for classical spin glass, reducing 
the initial many site problepi, to a one site problem. In order to obtain that, the partition function is expressed in 
terms of functional integrals^ using the anti-commuting Grassmann variables </>* s (t) and ^(t) (t is a complex time): 

Z = J D(4>*<t>) exp[A) + A pair + A SG + A f ie i d ] (5) 

where the action Aq in Eq. (||) is the free part of the total action, A pa i r is the pairing part and A$g is the spin part. 
Using time Fourier transformed quantities, the free part is given by 

Ao = J2 £ ( iu} + [$1 M&T M + til M^U M] ( 6 ) 

i u) 

and the A iield is 

Afield = ^J2J2 [tfrM&TM - tiiWuW] (7) 

i uj 

while the pairing part of the action in the static approximation becomes 

A P air^Af air = ^Y,Pi^PM ( g ) 



where /3 — — , T is the temperature and 



with Matsubara's frequencies u) = (2m + l)7r, m = 0, ±1, ... 

Indeed, functional integral tcchniqucs_within the static approximation have succeeded in studying phase transitions 
in conventional BCS superconductivity^ and in the presence of metal transition impuritiesEil. It has been possible 
to build a mean field theory in momentum space by using a particle-hole transformation. Thus, the problem became 
exactly resolved. 

The spin part of the action is written as 



where the spins variables are expressed in terms of Grassmann fields: 



(10) 



UJ 

sz{y) = Y t tii^ + v)^{u>) (12) 
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where v = 2irm in Eq. (p^). 

The Grand Canonical Potential can be obtained through the replica method by 

n = -lihn< z ">-- 1 . 

p n— »o n 



(13) 



The problem can be linearized and reduced to one site problem with the introduction of the auxiliary complex fields 
for the pairing part rj a (Q), rj*(Q) (r] a (0) = rj Ra (0) + irji a (0)) and for the spin partu E% a {v, v') = [B^ a {—v, -u')]* and 
Qapivt 11 ') — [Qa^( _y i where a and (3 are replica index and t = x, y, z. Thus the configurational averaged 



replicated partition function is 

(Z n ) ca = D( V * V ) J D(Q*Q) J D(R*R)exp- 



N 



EEi^-oi 2 + EEi^>X)i 2 

a/3 t<? a<f3 t 



8f3 2 J 2 
Nl3g 



EE \ Rt LM\ 

a t 

4 EM°)i 2+ln ( A ) 



(14) 



where K 



/ 2lt \ / 7T 



\NI3 2 J 2 J \Nf3g 
The function A in Eq. (Q is 



A = J W0)exp|A o + ^^[^(OK(O) + ^(OV(O)]+/3^E^(°) 



E E E O' v')sU-v)su-v) + E E E 

ioi(3 t<t' w' 



ia t 



ia.<j3 t vv' 



(15) 



Within the static approximation, together with the replica symmetric ansatz, the auxiliary fields of the spin glass 
part are giving by the following relationsEj: 



t^t' 

AI3 2 J 2 R Z 5 U . M ,8 U . Q 

RZ(u, is')=A[3 2 J 2 R5^,5^o 

4f3 2 J 2 Q z S u ^S ufi 



(16) 



where R, R z and Q z are now real parameters. 

We also take the spin glass order parameter transversal to the field Q — to work in a temperature region where 
both static approximation and replica symmetric ansatz are reliable. The sum over the replica index gives again 
quadratic forms which are also linearized by new auxiliary fields. Thus, the one site A function becomes 



A 



+ 00 



Dw 



Du x 



Du, 



+°° r _ _ 1 1 

Dv J D{4>*cj>) ex P L4 + A pair + p h ■ S]j 



(17) 



e~ v I 2 

where Dv = — dv and 



h ■ S = h z S z -r - 
h± = h x ± ih y 



h+S- + h-S-i 



(18) 
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The field h in the equation above is defined as: 



J 



v^/2{R z -Q z ) + w^2Q~ z + 



2J 



(19) 



h+ = JV2i?w 4 
h- = jV2Ru- 



The resulting problem obtained in Eq. ( |17| ) is equivalent to fermions in the presence of internal fields grj/2, grj*/2 
and h. The two first ones are the pairing fields related to the long range order where there is pairing formation in 
the sites. The later is a random gaussian field related to the replica the components of the spin part of the auxiliary 
fields introduced in Eq. (|l4|), namely, the non-diagonal Q z and diagonal R z which are parallel to H z , and diagonal R 
which is transversal to H z . 

In order to construct a mean field theory which can be solved, we introduce the following matrices: 



0J(-u) 



(20) 



(21) 



and 



G n (w) A 

/3/i+ G£(w) -A 

A* -g£(-"0 

-A* -GZU-w) 



(22) 



obtain 



In the equation above GjJ-(iv) — X) s =±(* w + Pt 1 + s /3frz) and A = — 77. Hence, for the A function in Eq. (|l7| ) we 



A 



+ OC 



Du,. 



Du[I(u;, tt x , tt„, u)] 



(23) 



where 



I(w,u x ,u y ,v) 



Dl 



) exp 



(24) 



The functional integral in Eq. (|24|) can now be performed. The sum over Matsubara's frequencies in the resulting 
expression should be done as in Ref. [l2|. The matrix formalism introduced in Eqs. (p0|), ( f2~lj ) and ( p^ ) produces a 
particle-hole transformation in the fermions of spin down. Therefore, the Grand Canonical Potential can be found 
from Eqs. © and © as 



o 1 / >+oc 

fi = 4/3 J 2 i? 2 + 2/3 J 2 R Z - 2(3J 2 Q Z + ^\ v \ 2 - - / Dwln^) 



(25) 



where 



uDu 



Dv 



cosh(/V) + cosh(/3|/i|) 



(26) 



and 



/3 M ' = v//3 M 2 + A 2 



(27) 
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with h already defined in Eq. (19). 

For Sp = sinh(/3|/i|)/(2/3|/i|) and 8 — h z /J, the functions R, Q z , R z and \r]\ (from now on we write rj instead of \rj\ 
can be determined from the saddle point equations that follow from Eq. (12a): 



+ OG 



Dw 



(28) 



2y/(R z -Q z ) J- 



Dw 



(29) 



Qz 



p 2 J 2 



Dw 



(30) 



(31) 



We have solved numerically the set of Eqs. fl2q), (j29|), ( |30D and (|31j) for the situation-where there is one fermion 
per site (in average) which corresponds to \x = due to the particle-hole transformation^. 

The validity range in temperature for this theory is the region where T > Tf implying that it is necessary to locate 
the Tf in this problem. 1C lis can be done by expanding the Grand Canonical Potential given in Eq. ([l4]) in powers of 
Q^p up to second orderlijO. The four spin correlation function that appears in the coefficient of the quadratic term 
can be related to the parameter R in the static approximation, where the averages are computed with an action where 
the auxiliary fields are given by Eq. ([l6]) and Q = 0. Thus the temperature Tf can be obtained from the condition 
that the coefficient of the quadratic term vanishes, so we get 



1 = i/3fJR 



(32) 



along the set of Eqs. (§|)-(§l|). 



III. PHASE DIAGRAM AND TRICRITICAL POINT 

The numerical solution of the equations (p8|)-([3l|) has allowed us to locate the transition temperature between the 
NP phase (rj = 0) and the PAIR phase (rj ^ 0) as a function of the pairing interaction strength g and the magnetic 
field H z (see Fig. (1)). The nature of the transition line has a complex dependence on both quantities. 

For instance, if H z =0 7 for high T and small g there is no long range order rj = (Q is always zero) that corresponds 
to the NP phase. On the other hand, for high T and high g the parameter 77 is nonzero. In this situation there is 
pair formation on the sites which corresponds to the PAIR phase. The transition between NP and PAIR phases is a 
second order type, which means that the order parameter rj goes continuously to zero as g decreases. That is shown in 
Fig. (2). However, if g is decreased bellow a particular value g = gt c (which corresponds to the temperature T — T tc ), 
the order parameter starts to display a discontinuous behavior indicating that the nature of the transition line has 
changed to a first order one. In the absence of the magnetic field the behavior of the other parameters are R=R Z and 
Q z = 0. 

If the magnetic field is turned on, the transition NP-PAIR only exists for larger values of the pairing strength g and 
the field tends to destroy the PAIR phase (Fig. (1)). The tricritical point (g tc , T tc ) is moved up and, as a consequence, 
the first order transition line exists over a larger interval of the pairing strength g. 

The {gtc, T tc ) point has been confirmed as a tricritical point by the expansion of the Gran Canonical Potential in 
powers of the order parameter 77 which defines the symmetry of the pairing phase: 



where 



pn = pn + A ir] 2 + A 2 i 1 i (33) 



Ai = / Dw^^ — -£ — -1 CUi 

J Q °° ue~%~du Dv[l + cosh(/3|/i|] 



G 



a 2 = [ + °° Dw 1 ~ I ue ^~ du J -™ Dv[1 + cosh ^ (35) 



The condition A\ = = A 2 (together with .the equations for R, R z and Q z ) gives the precise location of the 
tricritical point according to the known criteriaEx This location is shown in Fig. d3). Tricritical points have Jaeen 
already found in fermionic spin glasses for an Ising model with a pairing interaction in the half filling situationEj as 
well for an Ising model with charge fluctuationEj. 

The behavior of the parameters R, R z and Q z also changes in the presence of the magnetic field. Fig. (g) shows the 
results for H z = 0.00 where Q z — and R — R z . Figs. (||), (Q) and @ show the results, respectively, for H z /J/(8) 1 ^ 2 
= 0.25, 0.50 and 0.75 where one can see that the parameter Q z is no longer null while R^ R z . 



IV. CONCLUSIONS 



In this work we have investigated how the mechanism responsible for the BCS pairing formation in real space (PAIR 
phase) can be affected for a spin glass diagonal replica symmetric order parameter through a long ranged Heinsenberg 
model with a pairing interaction in presence of magnetic field in z direction. This model could be treated at mean 
field level within the static approximation and the replica symmetry ansatz, always with no transversal (to the field) 
spin glass ordering in order to remain in a region where the static approximation and replica symmetry are reliable. 
So, the initial problem has been reduced to one site problem with two different effective fields applied on, which are 
1] (rf ) and h. The former is a long range internal field related to a symmetry breaking which produces a pairing long 
range order and appears in the problem combined with the chemical potential to give an effective chemical potential. 
For the half filling case, just the pairing field survives. The later is a random gaussian field related to R, R z and Q z 
as defined in Sec. [n| As we can see in Eqs. (28) (pf|), the effective fields have to be solved coupled. That would be 
the mechanism through which the magnetic part of the problem always affects the PAIR phase, even above of the 
spin glass transition. 

By solving numerically the Eqs. (^8|)-(|3|), it has been possible to construct the phase diagram for several values 
of H z (Fig. (]l])). For high temperatures, it has been found a crossover from a continuous pair breaking symmetry 
to a sharp one where the PAIR order parameter r\ displays a discontinuous behavior as shown in Figs. (|| ^). The 
presence of the field H z creates an anisotropy which can be seen clearly in the behavior of the R, R z and Q z , shown 
in Figs. (||)-(H)- Therefore, the region in the phase diagram where there is a first order transition becomes larger. 
This behavior could be verified by positioning the tricritical point using the Eqs. (|34|) and (^5|). 

We remark that the our fermionic spin representation constitutes an important difference from Ref. |l5|. The 
consequences of different representations for the spins are clear if one compares our equations for the order parameters 
with those obtained in Ref. [l^. We also point out that the transition line between SG-PAIR is not attainable from 
the theory used in this work. 

Lastly, we have worked within the scope of the static approximation. The dynamic is absent in this formulation 
and certainly it is responsible for important effects at very low temperatures. Nevertheless, the static approximation 
is an upper bound of the theory from where the dynamic should be properly included. That would be subject for 
future work. 
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FIG. 1.: Phase diagram as a function of temperature and pairing coupling g/J for several values of H* where H* = 
H z / J/(8) 1 ^ 2 . Solid lines indicate second order transitions while dotted lines indicate first order transitions. Tricritical 
points are shown as filled squares. In the lower right corner it is shown the relation between H* and the NP-SG transition 
temperature Tf relative to Tt{H* = 0.00). The transition line between SG-PAIR is not attainable from the present 
approach. 
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FIG. 2. Dependence of the order parameters rj, R, R z and Q z for H* — 0.00 as a function of g/J for several values of 
temperature where T* = T/J. For high values of g/J the transition is clearly a continuous one as opposed to a discontinuous 
one for lower values of g/J. For H* — 0.00, Q z is always null. 
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